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The bulk parameters characterizing the energy of symmetric nuclear matter and the symmetry 
energy defined at normal nuclear density po provide important information on the equation of state 
(EOS) of isospin asymmetric nuclear matter. While significant progress has been made in determin- 
ing some lower order bulk characteristic parameters, such as the energy Eo{po) and incompressibility 
Ko of symmetric nuclear matter as well as the symmetry energy Esym{po) and its slope parameter 
L, yet the higher order bulk characteristic parameters are still poorly known. Here, we analyze the 
correlations between the lower and higher order bulk characteristic parameters within the framework 
of Skyrme Hartree-Fock energy density functional and then estimate the values of some higher order 
bulk characteristic parameters. In particular, we obtain Jo = —355 ± 95 MeV and lo = 1473 ± 680 
MeV for the third-order and fourth-order derivative parameters of symmetric nuclear matter at po 



and Ks 



-100 ±165 MeV, Js. 



: 224 ±385 MeV, h 



-1309 ± 2025 MeV for the curvature 



parameter, third-order and fourth-order derivative parameters of the symmetry energy at po, using 
the empirical constraints on Eo{po), Ko, Esym,{po), L, and the isoscalar and isovector nucleon ef- 
fective masses. Furthermore, our results indicate that the three parameters i?o(po), Ko, and Jo can 
reasonably characterize the EOS of symmetric nuclear matter up to 2po while the symmetry energy 
up to 2po can be well described by Esy,n{po), L, and Ksym- 
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I. INTRODUCTION 

The equation of state (EOS) of isospin asymmetric nu- 
clear matter, especially its isospin dependent part which 
is essentially characterized by the nuclear symmetry en- 
ergy, is important for understanding not only the struc- 
ture of radioactive nuclei, the reaction dynamics induced 
by rare isotopes, and the liquid-gas phase transition in 
asymmetric nuclear matter, but also many critical issues 
in astrophysics The nuclear matter EOS is conven- 

tionally defined as the binding energy per nucleon as a 
function of the density and a number of bulk parameters 
defined at normal nuclear density po are usually intro- 
duced to characterize the energy of symmetric nuclear 
matter and the nuclear symmetry energy. For example, 
the energy Eo{po) and incompressibility Kq of symmetric 
nuclear matter are the two lowest order bulk parameters 
for the EOS of symmetric nuclear matter while the sym- 
metry energy Esymipo) and its slope parameter L are the 
two lowest order bulk parameters of the nuclear symme- 
try energy. The bulk parameters defined at po provide 
important information on sub- and supra-saturation den- 
sity behaviors of the EOS of isospin asymmetric nuclear 
matter. 

While significant progress has been made in deter- 
mining some lower order bulk characteristic parame- 
ters of asymmetric nuclear matter, such as Eq{pq), Kq, 
Esym{po) and L [sl-fToj. yet the higher order bulk char- 
acteristic parameters are still poorly known. Actually, 
there is so far even not any direct experimental infor- 
mation on the third-order derivative parameter Jq of 
symmetric nuclear matter at po and the symmetry en- 
ergy curvature parameter Ksym- However, the higher 



order bulk characteristic parameters have been shown to 
be closely related to some important issues in nuclear 
physics and astrophysics, such as the determination of 
the isobaric incompressibility of asymmetric nuclear mat- 
ter [lllfl^ and the core-crust transition density and pres- 
sure in neutron stars Il3l4l5i. 



Theoretically, if the form of an energy density func- 
tional and its parameters are given, then the EOS of 
isospin asymmetric nuclear matter can be calculated and 
thus the bulk characteristic parameters at any orders 
can be easily obtained. Usually, the empirical values of 
some lower order bulk characteristic parameters, such as 
Eq{pq), Kq, Esymipo) and L, are used to constrain the 
parameters of an energy density functional. In such a 
way, while different energy density functionals usually 
predict similar results for the lower order bulk charac- 
teristic parameters, they could give very different predic- 
tions for the higher order bulk characteristic parameters. 
Therefore, the higher order bulk characteristic parame- 
ters can be sensitive to the energy density functional form 
and useful for constraining the energy density functional 
and its parameters. 

In the present work, we analyze the correlations be- 
tween the lower and higher order bulk characteristic pa- 
rameters of asymmetric nuclear matter within the frame- 
work of Skyrme Hartree-Fock energy density functional. 
Using the empirical constraints on the lower order bulk 
characteristic parameters and other macroscopic proper- 
ties of asymmetric nuclear matter, we then estimate the 
values of some higher order bulk characteristic parame- 
ters. 

The paper is organized as follows. We discuss the gen- 
eral properties of asymmetric nuclear matter in Section 
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im and then introduce the Skyrme Hartree-Fock energy 
density functional in Section Hill The resuhs and discus- 
sions are presented in Section IIVI A summary is then 
ffiven m Section El 



II. EQUATION OF STATE OF ASYMMETRIC 
NUCLEAR MATTER 



The EOS of isospin asymmetric nuclear matter, given 
by its binding energy per nucleon, can be expanded to 
2nd-ordcr in isospin asymmetry S as 



(1) 



where p = p„ + pp is the baryon density with p„ and 
Pp denoting the neutron and proton densities, respec- 
tively; 5 = {pn — Pp)/{pp + Pn) is the isospin asymmetry; 
^oip) — E{p,S — 0) is the binding energy per nucleon 
in symmetric nuclear matter, and the nuclear symmetry 
energy is expressed as 



-E-sym (p) 



1 d^Ejp.S) 
2! dS^ 



|<5=0- 



(2) 



The absence of odd-order terms in S in Eq. ([T} is due to 
the exchange symmetry between protons and neutrons in 
nuclear matter when one neglects the Coulomb interac- 
tion and assumes the charge symmetry of nuclear forces. 
Neglecting the contribution from higher-order terms in 
Eq. ([T|) leads to the well-known empirical parabolic 
law for the EOS of asymmetric nuclear matter, which 
has been verified by all many-body theories to date, at 
least for densities up to moderate values 0. As a good 
approximation, the density-dependent symmetry energy 
Esym{p) can thus be extracted from the parabolic ap- 
proximation of Esym{p) w E{p, (5 = 1) — E{p, S = 0). 

Around the nuclear matter saturation density po, the 
binding energy per nucleon in symmetric nuclear matter 
Eq{p) can be expanded, e.g., up to 4th-order in density, 
as 



Eoip)^Eoipo) 



2! 



Ji 



3! 



|x* + 0(x^), (3) 



where x is a dimensionless variable characterizing the 
deviations of the density from the saturation density po of 
symmetric nuclear matter and it is conventionally defined 
as 



Po 



3po 



(4) 



The first term Eq{pq) on the right-hand-side (r.h.s) of 
Eq. ([3]) is the binding energy per nucleon in symmet- 
ric nuclear matter at the saturation density po and the 



coefficients of other terms are 



Ko 
Jo 
lo 



9Po 



2d^Eo{p). 



dp^ 



lP=PO J 



dpS IP=P0' 



81p, 



jd-'Eoip)^ 
dp^ 



\P=Pa- 



(5) 
(6) 
(7) 



The coefRcient Kq is the incompressibility coefHcient of 
symmetric nuclear matter and it characterizes the curva- 
ture of Eo{p) at Po- The coefficients Jq and Iq correspond 
to the 3rd-ordcr and 4th-order derivative parameters of 
symmetric nuclear matter, respectively. 

Around the normal nuclear density po, the nuclear 
symmetry energy £^sym(p) can be similarly expanded, 
e.g., up to 4th-order in %, as 



Esymip) = Esymipo) + Lx + 



sym 2 
A 



+ - 



J. 



3! ^ 



4! 



-X 



2! 



(8) 



where L, Ksym, Jsym and /sym are the slope parameter, 
curvature parameter, curvature parameter, 3rd-order and 
4th-order derivative parameters of the nuclear symmetry 
energy at po, i-e.. 



L 



sym 



3po 
9Po 



dEsynijp) . 

dp ^"^"'^ 

2 <PEsyra{p) I 



9p2 



9p3 I 

^d^Esyn^p) I 



'sym 



81p, 



9p4 



\P=Po ■ 



sym, 



Jsyra and 



*sym 



(9) 
(10) 

(11) 
(12) 

characterize 



The coefficients L, 
the density dependence of the nuclear symmetry energy 
around the normal nuclear density po, and thus carry 
important information on the properties of nuclear sym- 
metry energy at both high and low densities. 

In the above Taylor expansions, we have kept all terms 
up to 4th-order in The 9 coefficients, namely, Eo{po), 
Ko, Jo, lo, Esymipo), L, Ksym, Jsym, ^sym, are theoret- 
ically well-defined, and they characterize the EOS of an 
asymmetric nuclear matter and its density dependence at 
the normal nuclear density po. Among these parameters, 
the lower order bulk characteristic parameters _Bo(po), 
Ko, Esymipo), and L have been extensively studied in 
the literature and significant progress has been made over 
past few decades. Based on the empirical constraints on 
the lower order bulk characteristic parameters and other 
macroscopic properties of asymmetric nuclear matter, we 
investigate in the following to what extend the higher or- 
der parameters Jq, Iq, Ksym, Jsym, and /sym can be con- 
strained within the framework of Skyrme Hartree-Fock 
energy density functional. 
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III. SKYRME-HARTREE-FOCK APPROACH 
AND MACROSCOPIC PROPERTIES OF 
ASYMMETRIC NUCLEAR MATTER 



In the standard Skyrme Hartree-Fock model, the nu- 
clear effective interaction is taken to have a zero-range, 
density- and momentum-dependent form (l6j . i.e., 



Fi2(R,r) = toil + xoPa)S{r) 

+ it3(l + a^3P.K(R)'5(r) 



+ t2{l + X2Pa)K' ■ S{r)K 

+ iWoiai+a^) ■ [K x S{r)K], (13) 

with r = ri — r2 and R = (ri -|- r2)/2. In the above, 
the relative momentum operators K = (Vi — V2)/2i and 
K' = — (Vi — V2)/2i act on the wave function on the 
right and left, respectively. The quantities P^r and ai de- 
note, respectively, the spin exchange operator and Pauli 
spin matrices. The a, to — t^, xo — x^ are the 9 Skyrme 
interaction parameters and Wo is the spin-orbit coupling 
constant. Within the standard form, the EOS of sym- 
metric nuclear matter can be written as 



Eoip) 



3^2 /37r 



10?71 



2\ 2/3 



3 ^ /3^2X 2/3 ^ ^ 

80 V 2 / 16 ' 



(14) 



with Qg = Ml + (5 -|-4a;2)t2- Furthermore, the symmetry 
energy can be obtained as 



Esym (p) 



1 f d^E ' 



6m 



24 



37r2 



2 \ 2/3 



1 



to(2a;o + l)p 



3^^ 



2\ 2/3 



— t3{2x3 + l)p' 



(15) 



ficient Gs, and symmetry-gradient coefficient Gy, i.e. 



to 
Xo 

h 
h 



Xi 



X2 



4a/(3po) 

3(2; - l)El;i^{po)/a - 1/2 
16/3/ [/9o^(7 + l)] 
-3y(7+l)i?^°?„(po)/(2/3)-l/2 

20c/ [9po(^-f)'] + 8Gs/3 
4(25G - 18D) _ 8(Gs + IGy ) 
9po(fcg)2 " 3 

12Gy - 4Gc - 



6L> 



/(3ti) 



(16) 
(17) 
(18) 
(19) 
(20) 

(21) 
(22) 



20Gv 

7-1 



4G5 



5(16G- 181?) 



■ipoiKY 



/(3i2) (23) 
(24) 



where the parameters G, D, a, /?, 7, £'s°^(po), and y are 
defined as 



G 



D = 



- rns.o ^0 



9 \ m 



m 



- 1 



-El 



kin 



Ko 



s,0 "'-u,0 / 

HG-^(i?°„-3i?o(po) 



lOG 



/? = ( 



Ko + ^Eo{po) 
E" 



EL 



4G 



EoiPo) - :;C) 



7 = 



kin 

3 3' 
Kq - 9Eo{po) + 5£;0i„ - 16G 
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^^L-4G 



-E's?m(Po) 



SEL 
E. 



sym(PO 
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'kin 



6G 



EsymiPo) 



D 



^ ' 3il-l)E%M 
withfc^ = (l.57r2p„)^/^^;0,„ = ^ ^3.' 



2D 



(25) 

(26) 

2G) 
(27) 

(28) 

(29) 
(30) 
(31) 



10?ri 



2/3 



2/3 



2/3 1 

Po ' and 



symvru; - (^Poj • Furthermore, Gg and Gy 
are respectively the gradient and symmetry-gradient co- 
efhcients in the interaction part of the binding energies 
for finite nuclei defined as 

i^grad = Gs(Vp)V(2p) - Gv [V(p„ - pp)]' /(2p). (32) 



with 9sym = 3^1X1 — ^2(4 -I- 5x2). 

As shown in Ref . Q , the 9 Skyrme interaction param- 
eters, i.e., a, to—t^, X0 — X3 can be expressed analytically 
in terms of 9 macroscopic quantities po, Eo{po), the in- 
compressibility Ko-, the isoscalar effective mass m* q, the 



isovector effective mass m 



s,Oj 

^ 0, Esymipo), L, gradient coef- 



IV. RESULTS AND DISCUSSIONS 

Based on the formulism in previous sections, one can 
now estimate the values of higher order bulk characteris- 
tic parameters, such as Jq and Ksym by analyzing their 
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FIG. 1: Jo from SHF with MSLO by varying individually L 
(a), Gv (b), Gs (c), £;o(po) (d), E,y^{p^) (e), (f), m^.o 
(g), ni;,o (h), po (i), and Wo (j). 



FIG. 3: Same as Figure [T] but for Ksyn 
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FIG. 4: Same as Figure [1] but for Jsym- 



FIG. 2: Same as Figure [T] but for /q. 



correlations with po, Eo{po), Kq, m* Q, ml Q, £'sym(po), 
L, Gs, and Gv within the framework of Skyrme Hartree- 
Fock energy density functional. As a reference for the 
correlation analyses below, we use the MSLO parame- 
ter set 0, which is obtained by using the following em- 
pirical values for the macroscopic quantities: po = 0.16 
fm-3, Eoipo) = -16 MeV, Kq = 230 MeV, m* □ = 0.8m, 
TO* = 0.7m, E^yn,{po) = 30 MeV, and L = 60 MeV, 
Gv = 5 MeV-fm^, and Gs = 132 MeV-fm^. And 
Wq = 133.3 MeV -fm'^ is used to fit the neutron P1/2—P3/2 
splitting in ^^O. It has been shown Q that the MSLO in- 
teraction can describe reasonably (the relative deviation 
from the data is less than 2%) the binding energies and 
charge rms radii for a number of closed-shell or semi- 
closed-shell nuclei. It should be pointed out that the 
MSLO is only used here as a reference for the correla- 
tion analyses below. Using other Skyrme interactions ob- 
tained from fitting measured binding energies and charge 
rms radii of finite nuclei as in usual Skyrme parametriza- 
tion will not change our conclusion. 

To reveal clearly the correlation of higher order bulk 
characteristic parameters with each macroscopic quan- 



tity, we vary one macroscopic quantity at a time while 
keeping all others at their default values in MSLO. Shown 
in Fig. [1] is the value of the higher order bulk character- 
istic parameter Jq. Within the uncertain ranges consid- 
ered here, the parameter Jq exhibits a very strong cor- 
relation with Ko and £'o(Po)- However, it depends only 
moderately on m* g and very weakly on po while it dis- 
plays no dependence on the other parameters £'sym(/5o)j 
L, Gs, Gv, iTT-lo, and Wq. From the empirical val- 
ues of So = -16 ± 1 MeV, Ko = 240 ± 20 MeV, and 
"^s — (0-8 ± 0.1)m, we can obtain an estimate of 
Jo = —355 ± 95 MeV. The results of a similar corre- 
lation analysis on the parameter Iq are shown in Fig. [21 
and it is seen that the parameter Jo also exhibits a very 
strong correlation with Kq and i?o(po) and the empirical 
values of £:o = -16 ± 1 McV and Kq = 240 ± 20 MeV 
lead to an extraction of Iq = 1473 ±680 MeV. Therefore, 
the higher order bulk characteristic parameters Jo and Iq 
can be determined within relative uncertainties of about 
27% and 46%, respectively, in the present analysis. 

In Fig. [3] we make a similar correlation analysis as 
in Fig. [T] for the symmetry energy curvature parame- 
ter Ksym- It is clearly seen that the parameter K^ym 
exhibits a very strong correlation with L and i?sym(po)- 
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FIG. 5: Same as Figure [T] but for Isyn 
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FIG. 6: (a) Energy per nucleon as a function of density for 
symmetric nuclear matter obtained by using Eq. ((3| including 
terms up to x'^i X^i ^'^d respectively, (b) Energy per 
nucleon as a function of density for symmetric nuclear matter 
obtained by using Eq. ((3} including terms up to with error 
bars due to the uncertainty of Jo. (c) Same as (b) but with 
error bars due to the uncertainty of /q. (i?o(po) = —16 MeV 
and A'o = 230 MeV have been assumed in the calculations.) 



However, it depends only moderately on m* q and m* q 
and very weakly on Eq{pq), Kq, and po while it is inde- 
pendent of other parameters G5, Gy, and Wq. From the 
empirical values of L 
30 ± 5MeV, 



60 ± 30 MeV and Esymipo) 



(0.8 ± 0.1)m, and m 



,•5,0 



(0.126 ±0.051)m [l3|, we can obtain Ksym = -100 ±165 
MeV. Furthermore, Figs. H] and [5] display the results 
from similar analysis as in Fig. [1] for the higher order 
bulk characteristic parameters Jsym and Igym, respec- 
tively. Similarly as in the case of Ksyrm they exhibit 
very strong correlation with L and £'sym(po)- The em- 
pirical values of L = 60±30 MeV, E, 
m* Q ~ (0.8 ±0.1)m, and m* n ~m. 
lead to the estimates of J, 



sym 



sym(po) =30±5MeV, 
„o = (0.126±0.051)to 
224 ± 385 MeV and 



I^y^ = -1309 ± 2025 MeV. These results indicate that 



the higher order bulk characteristic parameters of the 
symmetry energy, i.e., Kgym, Jsym and I sym J 

are largely 

uncertain based on our present knowledge in the stan- 
dard SHF energy density functional. 

Using the estimated values of the higher order bulk 
characteristic parameters, one can see how they influ- 
ence the EOS of asymmetric nuclear matter. Shown in 
Fig. ini is the case for symmetric nuclear matter. Fig. 




FIG. 7: (a) Symmetry energy as a function of density obtained 
by using Eq. ^ including terms up to x, X^i X^i a-nd x'', 
respectively, (b) Symmetry energy as a function of density 
obtained by using Eq. (|8j including terms up to x* with 
error bars due to the uncertainty of Ksym- (c) Same as (b) 
but with error bars due to the uncertainty of Jsym- (d) Same 
as (b) but with error bars due to the uncertainty of Isym-{ 
£'sym(po) = 30 MeV and L = 60 MeV have been assumed in 
the calculations.) 



[6t^a) shows the results obtained by using Eq. ([3]) includ- 
ing terms up to x^i X^-, and x^i respectively. One can 
see that Eq. ^ with terms up to [i-C-i the parabolic 
approximation] already gives a convergent prediction for 
the EOS of symmetric nuclear matter from about 0.5po 
to 1.5/90- The higher order terms of and x* with the 
characteristic parameters Jq and Iq in Eq. ([3]) signifi- 
cantly improves the approximation to the EOS at low 
densities and that up to about 2pQ. To describe reason- 
ably the EOS of symmetric nuclear matter above 2po, one 
needs to include higher order terms in x- These results 
arc consistent with those obtained in Ref. [l2| where the 
MDI interaction has been used. Figs. HKb) and (c) dis- 
play how the uncertainties of Jq and Iq affect the EOS 
of symmetric nuclear matter. It is seen that their un- 
certainties have only minor influence on the subsatura- 
tion density behaviors but significantly affect the EOS at 
higher densities (above about 2po). These results imply 
that the EOS of symmetric nuclear matter in densities up 
to about 2po can already be well described by the three 
bulk characteristic parameters Eq{pq), Kq, and Jq. 

Figure 7 displays how the higher order bulk character- 
istic parameters Ksym, Jsym and I sym influence the den- 
sity dependence of symmetry energy. Figure 7(a) shows 
the results obtained by using Eq. ([5]) including terms up 
to Xj X^) X^i and x^, respectively. It is seen that Eq. ^ 
with terms up to x^ can already give reasonably a con- 
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vergent result for the density dependence of symmetry 
energy up to about 2pQ. The higher order terms of 
and with the characteristic parameters Jsyvi and Isym 
in Eq. significantly improve the approximation to the 
symmetry energy up to about 3po- Figures 7(b), 7(c) and 
7(d) display how the uncertainties of Ksym, J sym and 
Isym affect the density dependence of symmetry energy. 
One can sec that the uncertainty of the parameter Kgym 
significantly affects both the sub- and supra-saturation 
density behaviors of the symmetry energy while the un- 
certainties of Jgym and Isym only have significant infiu- 
ence on the symmetry energy at higher densities (above 
about 2/5o)- These features indicate that the three bulk 
characteristic parameters £'sym(po), L, and Ksym essen- 
tially determine the symmetry energy with the density 
up to about 2pQ. 

V. SUMMARY 

We have analyzed the correlations between the lower 
and higher order bulk characteristic parameters of asym- 
metric nuclear matter within the framework of Skyrme 
Hartree-Fock energy density functional. Based on these 
correlations, we have estimated the values of some higher 
order bulk characteristic parameters. In particular, we 
have obtained Jq = -355 ± 95 MeV, /q = 1473 ± 680 
MeV, Ksym = -100 ± 165 MeV, Jsym = 224 ±385 MeV, 
and Isym = —1309 ± 2025 MeV using the empirical con- 
straints on Eq{pq), Kq, Esymipo), L, and the isoscalar 



and isovcctor nucleon effective masses. Our results indi- 
cate that the three bulk characteristic parameters Eq{pq), 
Kq, and Jq essentially determine the EOS of symmetric 
nuclear matter in densities up to about 2po while the 
three bulk characteristic parameters i?sym(po)j and 
Ksym well characterize the symmetry energy in densities 
up to about 2pQ. For higher density (above 2po) behav- 
iors of the EOS of asymmetric nuclear matter, the higher 
order bulk characteristic parameters become important. 

In the present work we have estimated the higher order 
bulk characteristic parameters only based on the stan- 
dard SHE energy density functional. It will be interesting 
to see how our results change if different energy- density 
functionals are used and these studies are in progress. 
The experimental information or model independent pre- 
diction on higher order bulk characteristic parameters of 
asymmetric nuclear matter is expected to put important 
constraints on the nuclear energy density functional form 
and its parameters. 
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